Crystal vibration limit of nanoscale dynamic etalon. Rough estimate 
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The paper analyzes the limits of nanoscale dynamic etalon caused by thermal vibration of crys- 
talline structure. It is shown that exactly sufficiently long one- dimensional defect of crystal lattice 
should confine the standard accuracy within angstroms. 
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One of the novel techniques of measuring objects of 
nanometer and subnanometer scales is related to creat- 
ing dynamic etalons based on the precise measurements 
the tunnel microscope probe movements []]. A linear 
dimension is transmitted from the etalon to a measur- 
ing system by means of a controlled movement of the 
etalon surface at a required distance and measurement 
of the above movement. In this case the real accuracy 
attained by the measuring system is limited by unpre- 
dictable vibrations of the crystal surface. The present 
paper is devoted to analysis of the fundamental limits 
of such system caused by thermal vibrations of crystal 
lattice. 

The thermal deviation of crystal atoms from the equi- 
librium states is analyzed within the adiabatic approxi- 
mation assuming the electron polarization to follow the 
relatively slow motion of atoms without delay. For molec- 
ular or covalent crystals like diamond, Ge, Si, GaAs the 
electron polarization effect is ignorable, whereas for ionic 
crystals (NaCl, CsCl, ZnS, etc.) this assumption is well 
justified for all the wave vectors of phonons except for 
a sufficiently small wave vectors k p ~ fl^/e/c (see, e.g. 
Ref. 0,13). Here fHs a characteristic frequency of opti- 
cal phonons, e is a crystal permittivity, and c is the light 
speed. In particular, for ZnS, where f2 ~ 6 x 10 13 s _1 
and e ~ 7, we get k p ~ 5 x 10 3 cm -1 . By contrast, ac- 
cording to the results to be obtained, actually phonons 
with wave vectors about t<l/a~2x 10 7 cm -1 (a is 
a characteristic lattice spacing) mainly contribute to the 
characteristic values of atom thermal fluctuations. 



(see, e.g., Ref. 0) 
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Here N is the total number of unit cells, the sum runs 
over all the possible wave vectors {k ^ 0} and all the 
modes, s = 1,2, ...,3a, describing different vibrations 
of atoms within one unit cell, m a is the mass of atom 
a, fi s (k) is the frequency of the vibration type labeled 

by the pair {k, a}. The eigenvectors u QS (k) = {uS(k)} 
(I = 1,2,3) of the corresponding vibration types meet 
the equalities 



u QS (k) • iw(k) = S SjS 



(2) 



where S s s i is the Kronecker symbol. Besides, the phonon 
frequency meets the equality O s (— k) = Q a (k), as well as 
the eigenvectors do, u QS (— k) = u as (k). 

In order to estimate the characteristic amplitude Sh of 
atom thermal fluctuations near the equilibrium states let 
us calculate the value 
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averaged over all the possible states of the equilibrium 
phonon ensemble with the temperature T and Hamilto- 



To estimate possible effects caused by the crystal sur- 
face states and its defects we analyze 1D-, 2D- and 3D- 
models. It should be noted that in the cases of ID- and 
2D-models the corresponding lattice is assumed to be 
embedded in the 3D-space, so its atoms actually perform 
3D-motion. 

Let the lattice unit cell contain a atoms and the value 
r na specify the displacement of atom a (a = 1, 2, . . . , er) 
in unit cell n = {n\, ri2, n^} (n$ € Z) from its equilibrium 
position R n = niai + ri2&2 + where {a^} are the 

crystalline spacing vectors (for the 2D-lattice U3 = and 
for the atom chain 77,2 — also). Then the corresponding 
operator f nQ and the conjugate creation-annihilation op- 
erators &k s , 6 ks of phonons are related by the expression 
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In this way we get 
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where the occupation number of phonons in the state 
Ik, under the thermodynamic equilibrium with tem- 
perature T (measured in energy units) is 



n s (k) 
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and M s (k) is some effective mean mass of atoms in the 
unit cell depending on the vibration mode that obeys the 
equality 



Due to equality J5J) the value M s (k) belongs to the in- 
terval M s (k) S [m m i n ,m max ], the boundary of the given 
interval are specified by the minimal and maximal masses 
of the unit cell atoms. 

In estimating the value of ((8h) 2 ) let us ignore the 
dependence of effective mass M s (k) on the wave vector 
k, i.e. set it equal to some constant, M s (k) = M e g- 
Within the adopted approximation sum (J5J is reduced to 
the integral over the possible values of the phonon energy 

e = mi-. 
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Here the density of the phonon states p(E) describes their 
distribution over the interval from zero to some maximal 
value £ max and is to be normalized to unity, 

"-Emax 

dEp(E) = 1 , 



(9) 



because the total number of states | k, s) must be equal 
to 3N a. Thereby estimate (0 can be represented as 
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where, as before, a denotes the characteristic value of the 
crystal lattice spacing. 

The upper boundary of the phonon energies -E max can 
be estimated using the characteristic frequency of the op- 
tical phonons f2 op t or what is the same the characteristic 
frequency of individual atom vibrations in crystal lat- 
tices, -E m ax ~ h$l opt (see, e.g., 0). For example, dealing 
with ZnS we have f2 op t ~ 6 x 10 13 s _1 , which corresponds 
to temperatures about Tp ~ M7 opt /fc,B ~ 460 K. 

To be specific in what follows we confine our consider- 
ation to the case of not too high temperatures, T < Tjj. 
Besides, setting M s g ~ 40 atomic units (for ZnS atomic 
masses of the compounds are Mz n 1=3 69 and Ms ~ 32) 
and a ~ 3A the following estimates of the cofactor in 
formula l(TT)|) 

h 2 
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is obtained. Therefore it is possible to ignore the effect 
of quantum fluctuations (the summand 1/2 in expres- 
sion O), so 
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If the phonon spectrum would be characterized by a 
single frequency f2, i.e. p(E) = S(E — Ml), then the 
desired estimate has taken the form 
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This expression actually corresponds to the model of a 
disconnected oscillator ensemble and the lower line of ex- 
pression (|13fl matches the estimate of the amplitude of 
classical oscillator affected by the white noise with in- 
tensity T. For T ~ ^,£1 and £1 ~ f2 op t the value of the 
given estimate is evaluated by expression (|1 1|) . It is ac- 
tually an estimate of the thermal vibration amplitude of 
a point defect in the bulk or surface of the crystal. So 
their effects on the accuracy of the nanoscale dynamical 
standard may be ignored. 

In order to get a more realistic estimate we will make 
us of the Debay model approximating the real phonon 
spectrum with solely acoustic phonons where f2 s (k) cx k. 
It is equivalent to approximating the phonon density p 
as follows 



p(E) = D 
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where D is the dimension of the crystal lattice under 
consideration, D = 1,2,3. Then, taking into account 
the adopted assumption about the value of temperature, 
T < To (recall that ksTjj := I? max := MId), we can set 
the upper boundary of integral (|12|l to infinity and, thus, 
rewrite expression (11 2t as 
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For the 3-dimensional lattice crystal integral (|15|l is re- 
duced to the following 
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i.e. again thermal vibrations of solid bulk has a minor 
effect on the accuracy of dynamical nanostandart. 

For the 2-dimensional (plane) lattice integral (|15|l has 
a formal logarithmic singularity that has to be cut off at 
small energies related to phonons with very long wave- 
lengths bounded by the crystal size or surface structures. 
It leads to appearing some cofactor L n in the expression 



(Shy 



(17) 



but its value cannot be too large, typically such effects are 
reduced to L n ~ 10. So, again, possible phonon modes 
localized at crystal surface which have be simulated with 
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the 2-dimensional lattice, has an insignificant effect of 
the standard accuracy. 

For the 1-dimensional chain the situation changes dra- 
matically. Integral H15|l has a singularity of the type 1/x 
and its long-wave vibrations give rise to significant devi- 
ation of atoms from the equilibrium positions, i.e. 




because the lower boundary of the phonon energy can be 
estimated as hc s /ell where c s ~ 10 5 cm/s is sound ve- 
locity in solids and I is the length of the chain. Setting 
ksT ~ Mlr> r*> hc s /a and I ~ 1 mm we get the conclu- 



sion that the latter cofactor in expression (|18|l gets the 
order of 10 7 . In fact, the considered oscillator chain simu- 
lating one-dimensional defects of crystal lattice is rather 
formal. Real atoms forming such defects embedded in 
solid bulk cannot deviate from their equilibrium positions 
for distances exceeding the lattice spacing. So we could 
expect that exactly sufficiently long low-dimensional de- 
fects limit the accuracy of the dynamical nanostandard 
on scale about angstroms on the side of thermal solid 
vibrations. 
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